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ABSTRACT: 

In  this  paper,  we  present  results  on  estimating  material  model  parameters  from  inverse  analysis  of  full-field 
deformation  data  obtained  with  a  prototype  of  a  novel  integrated  tool  consisting  of  a  digital  image  correlation 
system  and  software  for  data  analysis  and  parameter  estimation.  Such  a  tool  is  needed  for  characterizing  the 
properties  of  new  materials,  and  for  calibrating  and  validating  material  models.  The  stereo  microscope-based 
image  analysis  system  may  be  used  for  measurements  at  temperatures  up  to  750"C.  and  field  sizes  of 
approximately  1  mm.  The  Graphical  User  Interface  (GUI)-based  parameter  estimation  tool  integrates  image  data 
analysis  and  inverse  analysis  modules,  and  incorporates  features  for  interfacing  the  tool  with  commercial  finite 
element  (FEW)  packages.  Data  analysis  includes  filtering  or  correcting  for  rigid  body  translations  and  rotations, 
and  averaging  for  noise  reduction.  The  GUI.  together  with  a  micrograph  of  the  sample,  is  used  to  select  a  subset 
of  the  imaged  region  for  analysis,  and  for  specifying  sample  grain  boundaries  needed  for  developing  the  FEW 
model.  The  inverse  analysis  module  runs  the  FEM  simulation  under  the  experimental  loading  within  its  iterative 
loop,  using  the  downhil  simplex  method  for  parameter  estimation.  The  methodology  was  successfully  validated 
from  measurements  on  a  superalloy  sample. 


1  An  Integrated  Tool  for  Full-field  Measurements  and  Model  Calibration 

In  recent  years,  complex  multi-scale  material  models  have  increasingly  been  used  in  design  and  analysis  for 
various  engineering  applications.  The  scope  of  these  physics-based  models  is  varied,  dealing  with  stress-strain 
relationships,  defect  initiation  and  growth,  thermal  behavior,  fracture  mechanics,  etc.  The  purpose  of  attaining 
predictive  capabilities  using  the  material  models  is  to  estimate  the  performance,  safety  and  reliability  of  the 
components  that  are  developed  using  these  materials  and  operate  under  a  wide  range  of  conditions,  whether  it  is 
for  lifetime  prediction  of  existing  mechanical  components,  or  for  predicting  the  performance  limits  of  components 
designed  with  new.  engineered  materials.  However,  for  these  models  to  be  used  reliably,  they  need  to  be 
calibrated  and  validated. 

A  formal  overall  approach  was  presented  by  Johnson  el  a/.  [1]  with  the  goal  of  efficient  model  validation.  This 
approach  integrates  experiments  with  a  software  platform  where  optimization  techniques  use  FEW  models  of  the 
sample  and  the  data  to  estimate  model  parameters.  In  a  typical  experiment,  the  solid  sample  is  subjected  to  a 
sequence  of  step  changes  in  load  conditions.  A  series  of  images  of  a  sample’s  region-of-interest  (ROI)  is 
acquired  at  each  load  step.  The  images  are  then  precisely  aligned  and  evolution  of  the  displacements  is 
calculated  from  the  difference  of  each  of  the  images  relative  to  a  reference  image  following  extensive  image 
analysis  (1).  As  shown  in  Figure  1,  an  iterative  inverse  analysis  is  required  in  order  to  estimate  the  material 
parameters  of  the  model. 

We  have  developed  this  concept  into  a  novel  integrated  measurement,  analysis,  and  parameter  estimation 


: 


tool  that  enables  materials  engineers  and  researchers  to  make  deformation  measurements  of  materials  at  high 
magnification  and  at  high  temperatures,  then  analyze  the  data,  and  run  parameter  estimation  calculabons  using 
commercial  FEM  models  of  the  sample  in  iterative  loops  within  the  tool  In  this  section,  we  give  an  overview  of  the 
tool;  the  following  sections  detail  results  obtained  for  a  material  parameter  estimation  application 


We  have  designed  and  developed  a  new  high-magnification  stereo  microscope-based  Digital  Image 
Correlation  (DIC)  system  (Vic  3D)  that  can  measure  deformation  up  to  temperatures  of  750°C  with  sub-micron 
resolution.  DIC  is  a  non-contacting  measurement  technique  widely  used  in  experimental  mechanics  to  acquire 
full-field  deformation  measurements  [2H6).  The  technique  uses  a  random  dot  pattern  applied  to  the  test  article  to 
track  the  deformation  on  the  surface  (see  Figure  2).  During  the  experiment,  images  are  taken  at  intervals,  and 
digital  image  correlation  is  then  used  to  determine  the  deformation  in  the  image  plane.  In  order  to  relate  the 
displacements  found  in  the  image  plane  to  object  displacements,  a  calibration,  is  performed.  The  calibration 
permits  a  three-dimensional  reconstruction  of  the  object  shape,  as  well  as  the  determination  of  the  complete, 
three-dimensional  displacement  vector  field. 

Tima 
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Figure  2:  Tracking  of  objeci  deformation  using  Digital  Image  Correlation. 


The  underlying  concept  of  DIC  is  the  measurement  of  displacement  (and  possibly  displacement  gradients  as 
well  as  higher-order  terms)  at  a  point  in  an  image  by  considenng  the  gray  value  pattern  in  a  small  neighborhood. 
Given  a  second  image,  the  objective  of  the  method  is  to  find  the  displacements  (and  possibly  higher-order  terms) 
that  maximize  the  similarity  between  the  pattern  in  the  original  image  neighborhood  and  a  displaced  (and  possibly 
deformed)  version  of  the  same  pattern  in  the  second  image.  An  efficient  method  to  solve  this  problem  is  the 
iterabve  algorithm  proposed  originally  in  |3).  Let  I)  be  a  parameter  vector  that  contains  the  image  displacements 
and  possibly  higher-order  terms.  Let  /,  denote  the  gray  values  in  the  original  image  neighborhood,  and  g  (0) 


denote  the  corresponding  gray  values  in  the  image  after  deformation.  An  objective  function,  E.  is  defined  and 
minimized  using  standard  numerical  methods  to  obtain  the  best  parameter  estimate,  0. 


The  objective  function  could  be  based  on  the  sum  of  squared  differences  as  shown  above,  or  use  the 
normalized  cross-correlation  as  well  as  the  zero-mean  normalized  cross-correlation  criteria  for  improved 
performance.  The  resolution  capabilities  ol  the  DIC  method  have  been  shown  to  be  better  than  0.01  pixels  for 
displacements  with  negligible  bias  [5|.  |6J.  Since  the  digital  image  correlation  method  acquires  displacement  and 
strain  data  over  the  entire  specimen  area  being  imaged,  it  is  ideally  suited  to  measure  non-homogeneous 
deformation  fields. 
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Figure  3:  Functional  and  architectural  schematic  overview  of  the  tool  prototype.  The  images  1  and  2  show  the 
deformation  fields  from  FEM  computations  and  the  DIC  data,  respectively. 


Figure  3  shows  a  high-level  schematic  of  the  tool.  The  deformation  data  produced  by  Vic  3D  is  an  input  to  the 
Parameter  Estimation  Toolbox  (PET),  which  consists  of  eight  functional  modules  that  exchange  data  (mostly 
through  files)  as  shown  by  the  arrows  in  Figure  3: 

•  The  Model  Interlace  Module  implements  a  protocol  lor  exchanging  data  with  the  user-created  FEM 
model  (using  third-party  software)  of  the  sample  that  incorporates  the  material  model. 

•  The  Geometry  Delinltion  Module  allows  the  user  to  inleractively  identify  the  grain  boundaries  in  the 
form  of  polygons  using  a  digital  image  of  the  microstructure,  and  also  specify  regions  that  the  user 
wishes  to  exclude  from  the  parameter  estimation  computations,  e.g..  because  the  deformation  data  is 
not  sufficiently  accurate. 

•  The  Identifiability  Module  allows  the  user  to  run  the  FEM  model  and  test  the  sensitivity  of  the 
deformations  under  the  test  conditions  to  the  material  parameters. 

•  The  Data  Processing  Module  performs  extensive  pre-processing  of  the  data  from  Vic  3D  (including 
data  averaging,  transformations,  etc.)  to  allow  the  user  to  manage  the  large  amounts  of  deformation 
data  from  Vic  3D. 

•  The  Inverse  Analysis  Module  uses  optimization  algonthms  to  determine  the  optimal  values  of  the 
material  model  parameters  that  minimize  the  difference  between  the  measured  and  modeled 
displacements  while  running  the  FEM  model  in  an  iterative  loop. 

•  The  Results  Output  Module  incorporates  several  file  manipulation  and  plotting  routines  to  graphically 
display  results  of  data  analysis  and  parameter  estimation  computations. 

•  The  Graphical  User  Interlace  (GUI)  lAodule  serves  as  the  interface  to  all  the  other  modules,  and  it  ties 
the  various  modules  together  into  one  user-friendly  tool.  It  will  enable  the  user  to  develop,  control, 
and  examine  the  results  of  data  processing  and  the  parameter  estimation  process  on  a  Windows 
desktop. 


3 


•  The  Network  Module  allows  the  user  lo  control  the  process  using  the  GUI  Module  from  the  Windows 
machine  while  the  actual  computation  may  be  optionally  performed  on  a  different  high-performance 
computer  (Windows  or  Linux).  The  data  and  the  results  being  generated  are  available  to  the  user  for 
interaction  and  visualization  purposes  via  the  GUI.  The  modules  and  the  FEM  software  (shown 
above  the  dot-dashed  line  in  Figure  3)  can  run  on  a  separate  computer  while  the  remaining  modules 
run  on  the  host  Windows  computer. 

The  hardware-software  tool  is  innovative  in  many  ways.  The  high  magnification  stereo-microscope  system 
permits  deformation  measurement  with  sub-micron  resolution.  The  furnace  with  optical  access  enables 
measurements  at  high  temperatures  using  the  stereo  microscope.  The  parameter  estimation  software  has  been 
designed  to  be  modular  and  can  be  adapted  to  various  applications  relatively  easily.  For  example,  the  PET  can 
be  interfaced  with  a  wide  range  of  commercial  FEM  packages  and  various  data  formats.  It  can  incorporate  new 
optimization  techniques,  and  can  use  deformation  data  from  multiple  sources. 

2  Experiment  and  Data  Analysis 

As  a  demonstration  of  the  tool,  deformation  measurements  under  tensile  loading  were  performed  on  a  dogbone 
sample  made  of  single-crystal  nickel  superalloy.  The  dogbone  dimensions  were  50  mm  (length).  8  mm  (width.  3 
mm  in  the  neck  region)  and  1 .5  mm  (thickness).  The  rectangular  measurement  area  of  dimension  7.4  mm  by  2.8 
mm  was  located  at  the  center  of  the  sample  on  one  surface. 

During  the  experiment  chosen  for  validating  the  parameter  estimation  method.  388  image  pairs  were  acquired 
before  specimen  failure,  resulting  in  approximately  1.4  GB  of  image  data.  The  images  were  processed  using  Vic- 
3D.  and  approximately  60.000  data  points  per  load  step  were  analyzed.  The  Bunge-Euler  angles.  0;.  0  and  0j. 
were  used  to  characterize  the  orientation  of  the  crystal  axis  of  the  single-crystal  sample.  Sequential  rotations 
about  the  laboratory  coordinate  axes  (or  the  primary  axes  of  the  sample)  bring  the  former  into  alignment  v/ith  the 
material  axes  (or  the  principal  crystallographic  directions).  The  values  of  these  angles  were  measured  separately 
as  0,  =  8.6".  0  =  -3.2"  and  0>  =  -8.3". 

The  raw  DIC  data  was  averaged  at  each  FEM  mesh  point  using  an  exponential  weighting  function  that 
dropped  whose  shape  is  specified  by  the  user  of  the  tool.  For  this  sample,  the  displacement  data  at  each  of  the 
1071  grid  points  is  an  average  of  between  50-200  raw  data  points  (depending  on  the  measurement  locations)  for 
a  total  of  over  60.000  data  points.  The  data  was  then  filtered  to  remove  rigid  body  translation  and  rotation  using 
following  filter  at  each  load: 


F  = 


Here.  u0  and  v0  are  the  offsets  to  account  for  displacement  measurements  starting  at  a  non-zero  load.  The  third 
operation  is  the  rotational  transformation  which  corrects  the  rigid-body  rotation.  We  will  refer  to  the  filtered 
displacement  data  as  i/f.  v/.  and  w<.  The  filtered  v  profiles  are  shown  in  Figure  4. 

For  the  plastic  deformation  data  we  found  that  the  filtering  could  not  correct  the  distortions  in  the  y- 
displacements  because  these  distortions  are  not  caused  by  rigid  body  rotation,  but  due  to  the  large  localized 
deformation  in  u,  and  w;  near  the  top  right  corner.  We  need  to  use  the  v,  data  for  our  parameter  estimation.  Since 
the  geometry  and  loading  are  symmetric  in  the  FEM  model,  non-symmetric  deformations  as  seen  in  the  data  are 
not  predicted.  Hence,  we  chose  to  use  the  change  in  the  width  of  the  measurement  area  for  our  optimization  cost 
function  computations  instead  of  w.  This  width  change  data  was  still  too  noisy  for  use  in  parameter  estimation 
calculations.  To  reduce  the  noise,  we  computed  the  best-fit  slope  of  the  normalized  load.  F  vs.  scaled  average 
change  in  width  of  the  measurement  area.  Svo  .  which  are  defined  as: 


v  +  v0  -  mcan(  v) 


min 


M) 


[1 
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Here.  Ft  and  F?  are  the  loads  at  the  beginning  and  end  of  each  plastic  region.  The  scaled  change  in  width  at  any 
toad  with  the  offset,  ( Sv)r .  was  calculated  as: 


(*?), 


max 


The  scaling  ensures  that  the  slope  has  a  magnitude  of  the  order  of  unity  rather  than  the  very  large  values  that 
would  be  obtained  if  the  actual  loads  and  changes  in  widths  were  used. 
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Figure  4:  Effect  of  translation  and  rotation  on  load  vs.  elastic  v  data  at  three  points  along  the  vertical  centerline  of 
the  sample.  Left:  raw  data.  Center:  data  after  the  mean  value  of  v  has  been  subtracted.  Right:  final  filtered  data, 
vf.  after  performing  rigid  body  rotation  on  data  from  the  second  figure. 


3  Finite  Element  and  Material  Models 


3.1  Geometry,  Mesh,  and  Boundary  Conditions 

We  developed  a  3D  FEM  model  in  ADINA®,  a  commercial  sottware  product,  for  the  dogbone  sample.  The  mesh 
used  consisted  primarily  of  quad  bricks,  and  the  total  number  of  mesh  points  was  9894.  The  measurement  area 
was  meshed  separately  from  the  rest  of  the  dogbone  with  a  rectangular  grid  of  51*21  points  along  the  x  and  y 
axes.  This  separate  meshing  was  performed  for  convenience  to  obtain  a  higher  mesh  density  in  the  measurement 
area.  One  end  of  the  dogbone  was  fixed  (resulting  boundary  condition  fixed  u,  free  v  and  w)  while  the  other  face 
had  force  boundary  conditions  (equal  forces  at  all  the  nodes  on  that  face  that  add  up  to  the  total  applied  force, 
free  u).  Additionally,  v  and  w  were  fixed  on  the  end  points  of  the  longitudinal  centerline  on  one  surface. 


3.2  Elastic  Region 

The  single-crystal  superalloy  has  been  modeled  as  a  simply  orthotropic  material,  with  stresses  and  strains  in  the 
pnmary  (matenal)  coordinate  system  given  by  [7  : 


CU 


Here,  the  elasticity  stiffness  matrix.  C  is  a  6*6  matrix,  and  for  simply  orthotropic  materials,  only  6  out  of  the  36 
elements  are  independent:  E,(=E3) .  E3 ,  Gti .  G33(=G13) .  v,2 .  va(=v,j). 


3.3  Plastic  Region 

The  non-linear  plastic  region  was  modeled  using  Hill’s  model  for  yield  and  a  multilinear  hardening  rule  as 
described  in  this  section.  We  assume  that  the  sample  does  not  undergo  creep,  ie„  that  the  strains  are  only 
dependent  on  the  load  and  are  time-independent  at  a  given  load. 

The  Hill's  yield  model  has  six  parameters  which  are  the  yield  stresses  in  the  directions  of  the  matenal  axes  ( 1. 
2,  3)  and  the  yield  stresses  for  pure  shear  in  the  matenal  planes  (12,  23,  13),  i.e..  ar  t,  a*  3,  nf  3,  «ryi3.  ati3,  <r,  J3. 
respectively.  We  use  a  single  universal  plasticity  modulus.  for  modeling  hardening  |8).  Effective  Poisson 
ratios  along  the  matenal  axes  (denoted  by  subscripts  1 .  2.  and  3)  are  used,  and  these  are  related  via  a  Lankford 
coefficient  and  volume  conservation  as  shown  in  the  equations  below. 
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The  Lankford  coefficients  were  originally  defined  for  rolled  plates  [9].  Typically,  three  specimens  are  cut  from 
a  rolled  plate,  with  angles  of  0°.  45°.  and  90"  to  the  rolling  direction  (direction  a  ).  The  ratio  of  the  plastic  strain  in 
the  specimen  direction  transverse  to  a’  (direction  'b')  to  the  through-thickness  (direction  ’c'l  plastic  strain  is  the 
Lankford  coefficient.  In  ADINA’s  implementation  of  this  plasticity  model,  there  are  five  independent  parameters: 

nyi.ro,  r,s,  rut.  and 

The  elasticity  and  plasticity  models  used  would  treat  the  single-crystal  superalloy  as  a  bilinear  system.  The 
plastic  region,  however,  is  non-linear.  To  accommodate  the  non-linearity,  we  divided  the  non-linear  region  into 
three  piecewise-linear  regions,  as  shown  in  Figure  5.  which  plots  load  vs.  u  that  was  measured  on  the  single¬ 
crystal  sample  at  the  center  of  the  measurement  region.  The  three  locally-linear  regions  are  marked  Plasticl. 
Plastic  2  and  Plastic3.  with  end-points  corresponding  to  loads  of  3.68  kN.  3.72  kN.  and  3.75  kN.  respectively.  The 
deformation  continues  to  increase  with  small  increases  in  load,  and  additional  locally  linear  regions  can  be 
identified.  However,  for  purposes  of  demonstration  of  the  parameter  estimation  technique,  we  win  restrict 
ourselves  to  three  zones.  Hence,  a  set  of  Lankford  coefficients  is  needed  for  each  of  the  three  linear  segments 
that  approximate  the  plastic  region.  Of  course,  the  relationship  to  yield  stresses  via  the  Hill  parameters  is 
meaningful  only  for  the  first  linear  plastic  zone. 


X<dtspia:amfiflt&  u{mm|  rolatoe  laloft  odgo  ax=a.  xiQ1 


Figure  5  .  Left:  Load  vs.  u  data  at  the  center  of  the  measurement  region  for  the  single-crystal  sample.  Right.  Load 
vs.  v  data  at  x  =  0.  y  =  -1.35  mm  for  the  single-crystal  sample. 


4  Model  Parameter  Estimation  for  the  Elastic  Region 


The  generalized  cost  function,  J .  used  for  parameter  estimation  is  computed  as  the  two-norm  of  the  normalized 
error  between  the  model  prediction  and  the  expenmental  values  of  V,  averaged  over  the  total  number  of 
measurement  axes  (n)  and  FEM  models  ( m ).  Here  /  represents  either  the  (filtered)  measured  displacements,  or 
the  best-fit  slope  of  the  load  vs.  displacement  graph: 


./  = 


m  i-1 


1 - 

* 

% 

i 

K  .... 

The  subscripts  E  and  M  refer  to  the  expenmental  data  and  model  results,  respectively.  The  number  of  good  data 
points  (i.e.,  mesh  points  where  the  displacement  measurements  are  sufficiently  accurate)  in  the  f*  model  along 
the  /*  axis  is  N,r  Since  we  are  using  only  one  FEM  models,  m  =  1.  Also,  since  we  are  not  using  w  data,  n  =  2. 
Each  data  point  has  a  weight  WilA  that  can  have  any  value  between  0  and  1. 


After  performing  sensitivity  calculabons  for  parameter  identifiability.  it  was  found  that  only  twc  of  the  six 
parameters  in  this  orlhotropic  model  showed  sensitivity  to  the  measured  displacements,  and  hence  only  these  two 
(E,=  £Jf  E2)  could  be  estimated.  The  reason  for  this  limitation  is  two-fold.  First,  the  out-of-pane  w  data  was  too 
noisy  to  be  useful.  Second,  for  small  values  of  Bunge-Euler  angles.  n:  and  a<  are  non-zero  in  the  material 
coordinate  system  but  of  very  small  magnitude.  Consequently,  the  dependence  of  u  and  v  on  material 


6 


parameters  other  than  E?  and  vti (through  elements  of  Cj)  is  very  weak.  I(  the  material  axes  are  substantially  out 
of  alignment  with  the  global  coordinate  axes,  the  stress  vector  will  be  full,  and  thus  the  displacements  measured 
(or  equivalently,  the  cost  function  for  parameter  estimation)  would  have  significant  sensitivity  to  all  the  material 
parameters. 

Inverse  analysis  was  performed  using  the  downhill  simplex  method,  a  robust  search  algorithm  that  has  been 
successful  in  finding  the  global  minimum  of  various  non-linear  cost  functions  [10].  Table  1  shows  the  values  of 
the  Young  modulus  along  the  ( f.  2)  material  axes  starting  from  two  different  sets  of  trial  values  that  are  quite  far 
apart.  Figure  6  shows  the  load  vs.  displacement  {«/»)  graphs  from  the  FEM  model  using  these  optimal  parameters 
and  with  the  data.  Good  agreement  between  the  model  and  data  is  seen  in  the  elastic  stage.  For  the  other  four 
parameters,  we  use  nominal  values  for  this  material  of  n *  =  =  0.38.  im  =  0.3.  G,2  =  G.j  =  100  GPa.  G29  =  1 10 
GPa.  The  values  of  £,and  E^obtained  here  compares  well  with  the  effective  modulus  of  138.3  GPa  in  the  loading 
direction  calculated  from  the  gross  measurements.  On  a  Linux  workstation  powered  by  an  AMD  Opteron  248 
running  at  2.2  GHz  with  12  GB  of  RAM  machine,  the  optimization  computations  took  approximately  half  an  hour  to 
converge  (Figure  7).  Almost  all  of  the  time  is  used  for  the  cost  function  computations,  i.e..  running  the  ADINA' 
simulations. 

Table  1 :  Optimal  parameters  (elasticity  model)  for  the  single-crystal  sample. 


Run  1 

Run  2 

Initial 

Final 

Initial 

Final 

E,  (GPa) 

80.0 

129.5 

180.0 

129.5 

E,  (GPa) 

80.0 

141.6 

180.0 

141.6 

V"l,>  .to* 

Figure  6:  Load  vs.  V/at  four  points  along  centerline  (y= 0).  in  the  elastic  regime.  Circles  are  experimental  data  and 
the  straight  lines  are  FEM  model  results  for  optimal  material  parameters 


Figure  7:  Convergence  plots  for  the  estimated  material  parameters  of  the  single  crystal  sample,  and  the 
associated  cost  function.. 


5  Parameter  Estimation  for  the  Plastic  Region 

As  slated  earlier,  we  have  modeled  the  plastic  region  as  consisting  of  three  piece-wise  linear  segments  (marked 
Plasticl.  Plastic2  and  Plastic3  in  Figure  5).  each  with  its  own  universal  plasticity  modulus  and  Lankford 
coefficients.  We  first  estimate  the  parameters  for  Plasticl  zone.  The  elastic  properties  are  specified  as  the  values 
determined  from  the  elastic  analysis  of  Section  4.  In  order  to  implement  the  model  for  plastic  zones  2  and  3.  we 
assume  that  the  boundary  point  between  zone  1  and  2  is  the  yield  point  of  a  fictitious  material  with  elastic 
properties  Ea.  as  shown  using  dotted  lines  in  Figure  5.  Next,  we  estimate  the  elastic  properties  of  this  fictitious 
material,  and  the  plastic  properties  (yield  stresses,  universal  plasticity  modulus  and  Lankford  coefficients).  We 
discard  the  elastic  properties  and  the  yield  stresses  and  retain  the  Lankford  coefficients  and  the  universal 
plasticity  modulus  for  the  zone  Ptastic2.  This  procedure  is  then  repeated  for  the  third  plastic  zone. 

Sensitivity  analyses  showed  that  for  this  single  crystal  sample  under  the  imposed  loading  condition,  we  can 
estimate  only  three  of  the  five  parameters:  oy.j,  ro.  and  the  universal  plasticity  modulus.  £*V  The  other  two 
coefficients  are  assigned  the  following  values:  rti  =  0.2  and  rK  =  0.05.  These  approximate  values  were 
determined  from  the  measured  strain  ratios  in  the  three  directions.  The  plasticity  model  involves  both  the  yield 
stress  and  hardening  behavior  model  parameters  (Lankford  coefficients).  Consequently,  in  our  linearized  (slope- 
based)  approach  to  the  plastic  zone,  we  have  to  tie  down  the  plasticity  slope  to  the  yield  point.  To  this  end.  the 
cost  function  was  modified  to  incorporate  the  error  in  the  displacements  at  the  yield  load.  The  cost  function.  J. 
was  re-defined  for  plastic  zone  1  as: 
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Here  u.  is  the  filtered  x-displacement  at  yield  load.  This  term  was  weighted  by  0.25  to  reduce  the  potential  effect 
of  noise  in  the  u  data  on  the  optimization  accuracy.  The  choice  of  the  factor  of  0.25  was  reached  after  trying  a  few 
other  values  for  fastest  convergence.  V  is  the  best-fit  slope  to  the  load-displacement  plot.  The  subscripts  j  refers 
to  the  axis  (x  or  y).  k  represents  the  data  point  number,  and  E  and  <U  refer  to  experimental  data  and  model  results, 
respectively.  The  optimization  computations  took  approximately  12  hours  to  converge  for  Plastic  1  zone,  and  6-7 
hours  for  the  other  two  plastic  zones  on  the  AMD  Opteron  workstation  running  at  2.2  GHz  with  12  GB  of  RAM. 

The  optimal  parameters  for  the  three  plastic  zones  are  shown  in  Table  2  and  Table  3.  Figure  8  shows  the 
comparison  of  the  model  and  data  over  the  entire  load  range  over  which  parameter  estimation  had  been  carried 
out.  The  figure  on  the  left  shows  graphs  of  average  stress  vs.  strain  at  a  single  point  (center  of  the  sample).  The 
quality  of  the  model  fit  can  be  examined  more  closely  in  the  graphs  on  the  right  side  of  Figure  8  where  u,  has 
been  plotted  individually  for  the  elastic  and  the  three  plastic  zones.  The  fits  are  very  good  overall,  and  the 
optimization  convergence  is  excellent. 


Table  2:  Initial  and  final  parameter  values  for  the  first  plastic  zone. 


Run  1 

Run  2 

Initial 

Final 

Initial 

Final 

r„  (GPa) 

10.0 

13.7 

18.0 

13.7 

(Tvt(GPa) 

0.40 

0.77 

0.80 

0.77 

ro 

0.2 

0.016 

0.025 

Table  3:  Initial  and  final  parameter  values  for  the  second  and  third  plastic  zones. 


Plastic  Zone  2 

Plastic  Zone  3 

Run  1 

Run  2 

Run  1 

Run  2 

Initial 

Final 

Initial 

Final 

Initial 

Final 

Initial 

Final 

£“,(GPa) 

10.0 

6.4 

5.0 

6.4 

7.5 

3.7 

5.0 

3.7 

To 

0.05 

0.1 

0.25 

0.092 

0.1 

0.092 

Figure  8:  Left  Graphs  of  average  stress  vs.  strain  in  x  and  y  directions  at  various  x  points  along  centerline  (y=0). 
Circles  are  expenmental  data  and  the  straight  lines  are  FEM  model  results  for  optimal  material  parameters.  Right: 
ui  at  all  points  along  the  x  direction  for  elastic  and  plastic  regions  with  displacements  in  each  zone  shifted  to  pass 
through  same  starting  point. 
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6  Conclusion 

We  have  demonstrated  the  capability  of  a  new  hardware-software  tool  to  measure  full-field  deformation  data  and 
then  process  the  data  for  use  in  estimating  the  sample  s  material  model  parameters.  The  stereo  microscope- 
based  DIC  system  can  handle  field  sizes  of  about  1  mm  and  high  temperatures  up  to  750'C.  Measurements  were 
made  on  a  single-crystal  nickel  superalloy  sample.  The  parameter  estimation  toolbox  processed  the  vast  amounts 
of  DIC  data  that  was  generated,  and  interface  with  ADINA.  a  commercial  FEM  software  to  run  the  model  in  an 
iterative  loop  during  optimization  calculations.  The  elastic  and  plastic  model  parameters  converged  to  the  same 
values  which  are  typical  of  this  material.  It  was  noted  that  the  choice  of  sample  and  the  quality  of  the  data  can 
have  a  major  impact  on  the  parameter  estimation  results. 


Acknowledgement 

The  work  was  funded  by  Air  Force  SBIR  Contract  *?  FA8650-05-C-7251 .  The  authors  would  like  to  thank  Dr.  D.  A. 
Johnson.  Dr.  R.  John.  Dr.  M.  J.  Shepard.  2"3  LL  S.  M.  Visalli.  and  2nJ  Lt.  S.  Haenng  at  AFRL  for  technical 
discussions  and  for  program  support,  as  well  as  Prof.  M.  A.  Sutton  and  A.  P.  Reynolds  at  University  of  South 
Carolina.  Columbia.  SC  for  their  technical  suggestions. 


References 

[1]  D.  A.  Johnson.  W.  J.  Porter  III,  R.  John.  "Full-Field  Techniques  in  the  Study  of  Multi-Scale  Nonlinear 
Material  Behavior.”  Proc.  SEM  Annual  Conference.  Experimentation  and  Modeling  at  the  Nanoscate. 
Society  tor  Experimental  Mechanics.  Charlotte.  NC.  2-4  June  2003. 

[2]  T.  C.  Chu.  W.  F.  Ranson.  M.  A.  Sutton.  W.  H.  Peters.  ‘Application  of  Digital  Image  Correlation 
Techniques  to  Expenmental  Mechanics.”  Experimental  tAechanics.  25  (3).  pp  232  -  245  (1985). 

[3]  M.  A.  Sutton.  M.  Cheng.  W.  H.  Peters.  Y.  J.  Chao,  and  S.  R.  McNeill.  ‘Application  of  an  Optimized  Digital 
Correlation  Method  to  Planar  Deformation  Analysis."  Image  and  Vision  Computing.  4  (3).  pp  143-150. 
1986. 

[4]  J.  D.  Helm.  S.  R.  McNeill.  M.  A.  Sutton.  ‘Improved  3-D  Image  Correlation  for  Surface  Displacement 
Measurement."  Optical  Engineering  35:7.  pp.  1911-1920.  1996. 


" 


|5)  H.  W.  Schreier.  J.  R..  Braasch.  M.  A..  Sutton.  ‘Systematic  Errors  in  Digital  Image  Correlation  Caused  by 
Intensity  Interpolation Opt.  Eng..  39(11).  2000. 

[6]  H.  W.  Schreier.  M.  A..  Sutton.  ‘Systematic  Errors  in  Digital  Image  Correlation  Due  to  Undermatched 
Subset  Shape  Functions.*  Experimental  Mechanics.  43(2).  2002. 

[7]  R.  M.  Jones.  Mechanics  of  Composite  Materials.  2"’  Edition.  Taylor  and  Francis.  1999. 

[8]  ADINA  R&D.  ADINA  Theory  and  Modeling  Guide.  Volume  1:  ADINA  Structures  and  Solids,  p.  278. 
October  2005. 

[9]  W.  T.  Lankford.  S.  C.  Snyder.  J.  A.  Bausher.  Trans.  ASM.  42:1197.  1950. 

[10]  W.  H.  Press.  B.  P.  Flannery.  S.  A.  Teukolsky.  W.  T.  Vetterling.  Numerical  Recipes  in  C.  The  Art  of 
Scientific  Computing.  Cambridge  University  Press.  Second  Ed.,  1992. 


It) 


